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TECHNICAL MEMORANDUM X-64600 


AN ALGORITHM FOR SOLVING THE POPOV^ 
CRITERION APPLIED TO SAMPLED DATA SYSTEMS 

SUMMARY 


A systematic means is provided for analytically solving the Popov 
stability criterion that was developed by Tsypkin [1] for application to sampled 
data systems. It is developed in a manner that is amenable to solution by a 
digital computer, regardless of the order of the polynomials defining the 
numerator and denominator of the z-transform of the transfer fimction repre- 
senting the linear part of the plant. 

INTRODUCTION 


Lindorff [2] has presented a derivation of Popov’s stability criterion 
applied to sampled data systems which follows the development prepared by 
Tsypkin [1] . It provides a criterion for establishing sufficient conditions 
for absolute stability of a class of nonlinear sampled data systems described 
in Figure 1. The system must be capable of being separated into a single- 
valued nonlinearity represented by (pi - ) and a linear plant represented by 
G(s). If this system meets the conditions, 


(piO) = 0 , 

(11 

K > >0 for cr 0 , 

(21 


and G(s) represents a stable plant with at most one free integrator, then 
the T sypkin version of the Popov stability criterion is met if 


Re lG(zi 




(3) 


1. A portion of this report was presented as Paper No. G5 at the 
Third Annual Southeastern Symposium on System Theory at the Georgia 
Institute of Technology, Atlanta, Georgia, on April 5-6, 1971. 



NONLINEARITY 



Figure 1. Model of class of nonlinear sampled data systems. 


where z = e , T is the sampling period, and G(z) represents the 

z-transform of G{s) which is evaluated at z = i.e. , on the unit circle 

in the z-plane. 


As with continuous systems, the stability requirement on the linear 
part of the plant can be weakened to include open-loop unstable plants [3] . This 
is shown through a coordinate transformation. 


G, (sj 
tr 


G(sl 

1 + eG(s) 


(41 


where gain e is chosen so that the poles of the transformed system open-loop 
transfer G is) lie within the unit circle in the z-plane [2], This is equiva- 

i/i 

lent to requiring the nonlinear function cpia) to lie within the sector [e,k] 
for the original system to be absolutely stable. 

In determining the upper bound K of the sector, the real part of G(zl 
must be determined in some manner. This may be done by plotting G(z> in 
acomplex plane or by analytically determining the real part of G(zl and then 
finding its minimum value. Either procedure becomes tedious as the order of 
a system under investigation increases. 


DESCRIPTION OF PROPOSED ALGORITHM 


A systematic means for analytically solving the criterion of equation (31 
is provided, as a subsequent development to one suggested earlier by the 
author in a more primitive form [4] . It is developed in a manner that is 
particularly amenable to solution by a digital computer or electronic calculator, 
regardless of the order of the polynomials describing the nuiuerator and 


2 





denominator of G(z). It requires only that G(z) be expressed in terms of 
two polynomials in z ; i. e. , 


n 



The criterion of equation (3) may be recast in the form 

NiDi + (1 - 1 

+ (1 - > -- , 

where 


(5i 


( 6 ) 


n 

= Z T, X , 

j=0 •’ •' 

n 

Ni = Z r, Y . 

j=0 ■’ •' 

d 

Dl = Z 1, X, , 

j=0 •* ^ 

d 

^2 = Z ’ ( 7 ) 

j=o ■’ •' 

and 


/3 = cos wT . (8> 

The terms X. and Y. are developed by means of the recursive 
relationships, ^ ^ 


X. , - 2/3X. + X. , 
J+1 J J-1 


= 0 


5 


Y. , - 2j8Y. + Y. , 
J+1 3 J-1 


= 0 


(9a) 


3 



where 


Xo = Yi = 1 , 

Yo = 0 , 

and 

Xi = /3 , 

or by the alternate recursive relationships, 

Y. = X. + /5Y. . (9b) 

J+i J 1 

See Appendix A for the development of these relationships. 

Figure 2 represents a logic diagram demonstrating the ease with which 
a general program can be established to solve equation (6) for any plant meet- 
ing the conditions specified in the introduction. The linear portion of the plant 
may be of any order as long as it can be expressed in terms of a z-domain 
open-loop transfer function. See Appendix B for computer program. 

DERIVATION 

The derivation of equation (6) from equation (3) is brief. Defining z 
as 

iwT 

z = e (10) 

it follows that z ^ may be expressed as 

2^ = X^ + - /3^ . ( See Appendix A . ) (11) 

It is apparent from the definition of equation (10) that 
Xo = 1 , Yo = 0 , 

Xi = ^ , Yj = 1 . (12) 
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PRINT T, wT, 
-1/Re |g} 


NO 



Figure 2. Digital computer logic diagram. 
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One may use equations (12) and the recursion relations of equations (9) to 
generate values of and Y. for j > 1, Then, equation (11) maybe 

substituted into equation (5), letting z = z, 


G(z) 


iwT 

z=e 



( 13 ) 


Taking the real part of equation (131 yields the left half side of equation (6) , 


EXAMPLE 


Several examples show how the proposed algorithm is implemented. 
Consider an integral control servo system [5] shown in Figure 3. The cor- 
responding z-transform of the linear portion of the plant may be put in the 
form of equation (5) and is described in Table 1. By using Table 1 as an input 
to the digital program, the left side of inequality (6) may be solved as coT 
varies in value from zero to tt for selected values of T. The results are 
printed out, and the minimum value of 

- 1/Re / G(z) 



is selected as the maximum value which K may possess. The results, repro- 
duced in Table 2, agree with those of Reference 5. The same algorithm has 
been used to find stability limits on K for two other examples of Reference 5 
(Bertram’s example and the integral-control example with digital compensa- 
tion added) . The results have been identical. 



Figure 3. Model of integral-control servo system. 
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TABLE 1. TRANSFER FUNCTION COEFFICIENTS 
FOR INTEGRAL-CONTROL SERVO SYSTEM 


j 



0 

1 - - Te-^ 

-T 

e 

1 

-T 

T - 1 + e 

-T 

-(1 + e ) 

2 

0 

1 


TABLE 2. STABILITY LIMIT ON K 



Sampling Period T 

1 

2 

3 

4 

5 

K< 

0.666 

0. 5 

0.4 

0.333 

0.286 


CONCLUSION 


A general computer program can be developed to determine the suffi- 
cient condition for stability of a class of nonlinear sampled data systems in 
terms of a sector within which the nonlinearity must lie. The criterion used 
in developing the computer algorithm is Tsypkin's development of the Popov 
stability criterion. Work is now in progress to investigate a similar algo- 
rithm for Jury's and Lee's [5] less restrictive stability criterion, 


Re{G(zl[l + q(z - 1)]} + 


K' |q! 


K 


(z - l)G(z) 


^ 0 


(14) 


where 


-00 :< 


dcj){cr) 

da 


< K' 


5 


and q is a real number. 
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APPENDIX 

DERIVATION OF RECURSION RELATIONSHIP 


The two sets of recursive relations, equations (9a) and equations (9b), 
are developed from the definition of z, 

iwT ^ ^ 

z = e - cos coT + 1 sin wT . (10) 

Letting 

/3 = cos coT , (8) 

z may be rewritten as 

z = /3 + i\li - . (A-1) 

If z is raised to the power, it may be written as 

= cos jcoT + i sin jcoT . (A-2) 

Two identities [6] for cos jcot and sin jcot may be rewritten in terms of 


io-3) 


cos jcoT = 


I '-'CT’XTiT)'' 


-2p-3^j-2p-2 

(A-3) 


and 




smjc.T=Vrr^ 2 . (A-4) 

p=o V P / 


where p is an integer. If sin jcoT and cos jcoT are defined as 


2. I ani indebted to Mr. Hans Hosenthien, Chief, E&D Analysis Office, 
Astrionics Laboratory, NASA, Marshall Space Flight Center, for the develop- 
ment of the alternate set of recursive relations of equations (9b), 
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cos jcoT = X. 


J 


and 


sin jwT = Y.n/TT , 


(A-5) 


examination of equations (A-3) and (A-4) shows that and are poly- 
nomials in )S. K equations (A-5) are substiiuted into equation (A-2), the 
resulting expression becomes of the form 


= X. + iY.N/i - . 

J } 


(11) 


The recursive relations of equations (9a) or (9b) may be derived by 


... -j+1 

writing z as 


-j+1 -j- 

z = z z 


= ^X.'+ iY.N/l - + is/l - ^2 ^ 

= X.i8 - Y.(l - + Y./3 ^n/TT^ 


= X. , + iY. . , 

J+1 J+i 


yielding the recursive ’ tions 

X. = X.i3 - Y.(l - /32) 
J+1 J J 


and 


Y. , = X. + Y.i3 . 
J+1 J J 


From the definition of z, 

Xo = = 1 , Yo = 0 , Xi = /? . 

Using equations (9b), one may readily derive equations (9a). 


(9b) 


(A-7) 
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APPEND l-XB 

DIGITAL COMPUTER PROGRAM DESCRIPTION 

The Popov-Tsyidtin digital program was written in Fortran IV, ai^d com- 
puter runs were made on an SEl^ 840 digital computer. Using the Tsypkih ver- 
sion of Popov's criterion, the program calculates and prints the negative 
reciprocal of the real part of the z-domain transfer function representing the 
linear portion of the plant. Values of the negative reciprocal are calculated 
for incremental values of the independent argument coT and scanned; the 
smallest value is selected as the upper bound of the nonlinear (or Popov) 
sector. The following mathematical model is described^by equations (3) and 
(6) through (9a), and is the basis for all program calculations (see program symbol 
table for relations of mathematical symbols and program symbols). 


Program Symbol Table 


Math Symbol 

Program Symbol 

Description 

Xo 

X(l) 

Initial value for X. 

Yo 

Y(l) 

Initial value for Y. 

Xi 

X'2) 

Initial condition for X. 

Yi 

Y(2) 

Initial condition for Y. 


AIN 

Increment denominator which is divided 
into RA to form the omega increment 
value. 


RA 

Range value of omega. 


WI 

Initial omega. 

T 

T 

T parameter which should be set to 1. 0 
and associated with wT limit values. 


DW 

Delta omega which value is used if 


omega equals 0. 
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Math Symbol 

i Program Symbol 

Description 

d 

ID 

Upper limit on Dj and D 2 summation and 
has value d + 1, 

n 

N 

Upper limit on Nj and N 2 summation and 
has value n + 1, 

K' 

KP 

Program flag which when set to 0 causes 
a Popov-Tsypkin run to be made. This 
relates to a program modification now 
under development to solve equation (14). 


See program inputs for all other values. 


Program inputs 

Card 1 — Reads X(l), Y(l), andY(2) under a 4F15. 8 format. These 
are the initial X and Y values. 

Card 2 — Reads AIN, RA, and WI which are the increment value of 

omega, the range of omega, and the initial omega. Reads under 
a 4F15.8 format. 

Card 3 — Reads T under an FI 5. 8 format. 

Card 4 — Reads DW, the delta omega value, under an Fi5.8 
format. 

Card 5 — Reads ID, the value for d , under an 15 format. 

Card 6 — Reads N, the value for n, under an 15 format. 

Card 7 — Reads KP, the value which when set to 0 causes the program 
to make a Popov-Tsypkin run using the initial equations to 
compute Re {G(Z)} and its negative reciprocal. When 
KP = 1, the program omits this phase and computes a K' 
equation set; where KP = 2, the above options are omitted 
and the program solves a quadratic K equation. For the 
Popov-Tsypkin option always set KP = 0, 
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Cards 8 through 11 contain the inputs involving options other than the Popov- 
Tsypkin equation solutions that were part of the program study and are included 
in the present version of the program. These inputs are included here only 
for the programmers information. 

Card 8 — Reads IQ, the number of Q elements, under an 15 format. 

Card 9 — Reads IK, the number of K elements, under an 15 format. 

Card 10 — Reads array Q, containing all values of q, under a 
4F15. 8 format. 

Card 11 — Reads array AK, containing all K elements, under a 4F15,8 
format. 


Sample Run Case 

The following sample case was run successfully on an SEL 840 digital 
computer. The appropriate gamma and eta functions listed in Table 1 of the 
computer printout were used in the Popov-Tsypkin program with the input 
data listed here and described under the card input data. The input data will 
always be followed by the printed variable values for T, coT, Re {g(Z)}, 
and K < -l/Re (G(Z)}. 
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Q::- e » 5 » # * !! * » » '^0P3V " rSV^^KI^ P-S0GRAH <> 9 ff » * 9 9 9 9 9 1 

dimension O'i^O) . AK(20) 2 

UrMf-NSI0N Xi 20 ) . Y ( 2u ) ,G(20 ) .E(20 ) 3 

OlMhNSlai'J IAl.A(20> 4 

TRACt 9999 5 

«E aO ( 4. 10 ) X( 1 ) » Y ( 1 j , Y< 2 ) 6 

READ (4.1U) AIN.RA.RI 7 

read (4.11) I 8 

READ (4,11) UW 9 

reap (4,12) id 10 

READ (4,12) N 11 

read (4.12) KP 12 

read (4,12) IQ 13 

READ (4,12) IK 14 

read (4,10) (Q< 1 ) , 1=1, pj) 15 

READ (4, 10) ( AK( I ) , 1 =1, IK) 16 

READ (4,120i IAEA 17 

10 F0RMAT (4riP,8) 18 

11 F0RHAT (115.0) 19 

,12 i-0RMAr (15) 20 

120 Format (20AA) 21 

WRITE (5,12U) IAEA 22 

WRITE (5,20T 23 

20 F0RMAT (5X, IEHINPUT DATA//) 24 

WRI 'E (5.21) X(i),Y(l),Y(2) 25 

21 Format (lOX./HXd) = , F l O . 5 . 3X. /HY ( 1 ) = ,F1U.5,3X, ;hy( 2) = ,F10.t>) 26 

WRITE (5,26) AIW.RA.WI 27 

26 Format ( lOX . ISH^MEGA increment = ,f 10 .5. 3x, 14H0M6GA RANGE = ,f'l0.5 28 

l,.iX,16H0MEGA INITIAE = .F10.5) 29 

write (5.22) ID 30 

22 F0RMAT (1UX.4HD = ,1R) 31 

WR I i E ( 5, 23> N 32 

23 F0RMAT (10X.4HN = ,15) 33 

WRpE (5,24) T 34 

24 F0RMAT (10X.4HT = ,H5.8) 35 

WRITE (5,27) DW 36 

2/ F0RHAT (lOX, 14M0EET A 0MtGA = ,F15.6) 37 

WRITE (5,50;i) KR.IQ.IK 38 

500 F0R.YAT (10X.5HKP = ,I5,5X,5HTQ = ,15.5X,5HIK = ,15) 39 

Write (5,50i) ( i , q( i ) , i = i. k^) 40 

501 F0RMAT (10X,2HQ(, 12.4H) = ,FiO,5) 41 

WRITE ( 5.502 ) ( I , AK( I ) , 1=1, TK) 42 

502 Format (10X,2HK(, !2,4M) = ,H0.5) 43 

W=0,0 44 

C«99F0RM gamma and eta function 45 

El = bXP(-’r> 46 

G(l,) = (1.0-ei)-T»El 47 

G(2)=T*(1,0-E1) 48 

E(l)=El 49 

£( 2 ) =-( 1. 0*tl) 50 

6(3) =1.0 51 

C9-9«F0Rm REAE R00I FU-iCTIIAiM 52 

IN=AIN 53 

IC=IN+1 54 

VA=»A/A1N 55 

ir(KP» 1)4010,4 01, 401 56 

4010 C0NTINU6 57 
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WRITE (5,41) 58 

41 F0RMAT (IIX. 1HT.18X.2HWT , 14'>(,8HRb(G{2) > , 4X.;^3HK LESS THAN -1/RE(G( 59 

li))) 60 

401 C0NTlNUb 61 

IF (KP-1) 39,301,390 62 

39 U^NTlNUfc 63 

D0 40 1=1, lU 64 

A I = i 65 

W = ,v'l+VA<fAI 66 

IF (W.EQ.O.U) W=DW 67 

WT=WoT 68 

8 = C'3S(Wn 69 

X(2)»B 70 

IF ( ID.LT.3) G0 T0 41C 71 

D0 30 J=3, lU 72 

X(j)=2,U»8*X(J-l)-X(J-2) 73 

Y( J)=2.0»3*YIJ-1)-Y( j”2) 74 

30 C0nUNU& 75 

call sum (N. id, IP,X. Y,E,G, AN l, AN2,D1,D2) 76 

G0 T0 300 77 

410 C0NriNUt 78 

Call sum (N, id, IP.X. Y,E-G,AN1,AN2,D1,D2) 79 

300 C0NTINUE 80 

ANUM = AN1»01+ (1 . 0-8*»2 )*AN2*D2 81 

DFm = D1«D1+(1. 0-8»*2)«02-"'D2 82 

RE=ANUM/DEN 83 

REC=-1.0/R£ 84 

write (5,31) T,WT,R£,REC 85 

31 F0RMAT (3X,E15.8.5X.£15.8,5X,E15,8,5X,E15.8) 86 

40 C0NTINUE 87 

G0 10 4200 88 

301 C0NT1NUE 89 

D0 380 I J=l, IK 90 

00 380 11=1,10 91 

D0 380 1 =1. IC 92 

AI=1 93 

Wsl|I+VA*AI 94 

IF (W.EQ.0.0) W=DW 95 

WT=W»T 96 

8 = C3S( W r ) 97 

X(2)’'B 98 

8M1=1.0-B 99 

8S0=?1 . O'B^B 100 

IF (ID.LT.3) G0 T0 4i00 101 

D0 3000 J=3, ID 102 

X( j)=2,0»8*X(U-l)-X( J-2) 103 

Y< J)=2.0«8*Y( J*l)-Y( J”2) 104 

3000 C0NTINUE 105 

call sum (N, id, IP.X . Y ,E,G, AN l, AN2,D1,D2) 106 

G0 T0 3111 107 

4100 C0NTINUE 108 

call sum (N, 10, IP.X. T.E.G, ANl, AN2,Dl.i)2) 109 

3111 C0NT1NUE 110 

RK=1.0/AK(1j) 111 

D8N = ANl»Dl-*'k5SQ#AN2»D2 112 

08=1. 0-0( II )»8M1 113 

8Q = Q(1I)»8SQ U4 

DN=AN1«D2-AN2*D1 115 
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DB = U1«D1 + 8SU»D2*D2 U6 

T0P = DBN»Q8-!'8O«D!'J*RK*D8 U7 

Aa8=Qnn'*8Mi us 

8n = AN1»AN1-»BS0»AN2»AN2 U9 

B0T = AOB^BiM 120 

AKP=T0P/80T 121 

WRITE IJ» AK( IJ> , I I ,0( |I ),WT» AKP 122 

383 F0RMAT (3X,2HK( . I2»4H) = , £15 , 8 , 2X » 2HQ ( , 1 2 . 4H ) = , El5 . 8 . 2X , 5HWT = 123 

l,Elb.8,2X, 7HK PR = ,£15. 8) 124 

380 Continue 125 

G0 T0 4200 126 

390 C0NTINUE 127 

D0 399 11=1. IQ 128 

00 399 1=1,10 129 

AI=i 130 

W = wUVA«AI 131 

IF CW.EQ,0,Q) W=DW 132 

WT=W*T 133 

8 = C0S(Wn 134 

X<?)=8 135 

8M1=1,0-B 136 

8SQ=1.0-B«B 137 

IF UD,LT,3) G0 T0 4101 138 

D0 3001 J=3. ID 139 

X(,J)=2.0«B»X( J-D-X( J-2) 140 

Y( J)52.0'->8*YCJ-l)-YU-2) 141 

3001 Continue 142 

CALL SUM (N, ID, IP,X, Y,e.G,ANl,AN2,Dl,02) 143 

G0 T0 3112 144 

4101 C0NJTIMUE 145 

Call sum CN. id. IP,X, Y,E,G,AN l,AN2,Dl,n2) 146 

3112 C0NT1NUE 147 

AAsQC I I )*8M1»( AMl»ANi + BSQ»AN2*AN2) 14 8 

B8=-( ANl»Dl+8SQ»AN2»D2)»(l,0-on I 149 
CC=-(D1»D1'*BSQ*D2«D2) 150 

DISC=S0RT (Bb«88-4, 0«AA*CC) 151 

AK1=(-BB+DISC)/(2.0*AA) 152 

AK2=(-B8-DISC)/(2,0«AA> 153 

WRIIE (5,3990) ( I I ,Q( 1 1 ) , AKl, AK2, WT ) 154 

3990 F0RMAT C3X,2HQ( , I2,4H> = , E l5 , 8, 5x , 5HK 1 = , E15 , 8', 5X , 5HK2 = ,E15.8, 155 

15X.5HWT = ,hl5,8) 156 

399 C0NiTlNUE 157 

4200 C0NflNUb 158 

WRITE (5,42) 159 

42 F0RMAT (25X‘. 14HCASE C0MPLETED) 160 

WRITE (5,43) 161 

43 F0RMAT (IHD 162 

Q0 T0 1 163 

EMD 164 


PATH 

00000 Q 

00050 AK 

00120 X 

00170 Y 

00240 G 

00310 E 

00360 lALA 

UNRF S9999 
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(} 0 4 0 

SI 

0 3323 

T$1003 

00t576 

$10 

suap 

SUM 

.!3251 

AIN 

0 3 524 

IP 

032^53 

RA. 

0 3325 

ani 

0 32b5 

W1 

0332 7 

AN2 

00601 

$11 

03331 

U1 

0 3? 57 

T 

0 3333 

U2 

0 3261 

DW 

01556 

$300 

0 0 6 0 4 

$12 

03335 

ANUM 

0 3?63 

ID 

0333 7 

DEN 

03264 

N 

03341 

RE 

03265 

KR 

0 3343 

RfcC 

03266 

13 

01661 

$31 

03267 

IK 

Q3217 

$4200 

03270 

I 

02457 

$380 

03271 

TSIOOO 

03 3 4 7 

I J 

00606 

$120 

03350 

II 

0 0 624 

$20 

0 3351 

8M1 

D0647 

$21 

03353 

8SQ 

00 704 

$26 

0 213 7 

$4100 

Q0 743 

$22 

02112 

$300 0 

00 757 

$23 

02156 

$3111 

00 7 /3 

$24 

03355 

RK 

01010 

$27 

03357 

DBN 

01036 

$500 

03361 

QB 

0110 0 

$501 

03363 

BQ 

01136 

$502 

03365 

ON 

03272 

W 

0 3367 

DB 

032/6 

Bl 

03371 

T0P 

S‘IBR 

EXP 

03373 

TS2001 

03 500 

T$2000 

0 33 75 

AQB 

03304 

IN 

033 77 

8N 

03305 

IC 

0 3401 

80T 

03306 

VA 

1)34 03 

akp 

01266 

$4010 

02430 

$383 

U1316 

$401 

03202 

$399 

012 7 5 

$41 

02/32 

$ 4 1 0 1 

01330 

$39 

02/05 

$3001 

01/05 

$301 

02/51 

$3112 

02502 

$390 

03405 

aa 

016/4 
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-.28302418E 

00 

-.2/896807E 

00 

-.2/496443E 

00 

■..271012686 

00 

".26711216E 

00 

-.26326229E 

00 

-.25946245E 

00 

-..25571208E 

00 

■»,2520105/E 

00 

-.24835733E 

00 

-.24475181E 

00 

•..24119344E 

00 

•..23768165E 

00 

-.23421589E 

00 

-.23079561E 

00 

•..22742027E 

00 

-.22408935E 

00 

-.22080231E 

00 

-.2.L755862E 

00 

■..21435779E 

00 

-,2lll993')E 

00 

•>.2O0O8264E 

00 

-.2D500734E 

00 

-.20197289E 

00 

19897881E 

00 

» . 19602464E 

00 

■..1931099UE 

00 

” . 19023413E 

00 

-.18739688E 

00 

-.18459769E 

00 

-.181836136 

00 

-.179111756 

00 

-.1.'642412E 

00 

-.173772836 

00 

-.1711574'->E 

00 

- . 16857 7566 

00 

1660327 /E 

00 

-.1635226/6 

00 

« , 1 D 1 0 4 6 8 ■' 6 

0 0 

- . 158604986 

00 

■». 156196616 

00 

-.153821386 

00 

- . 15147894E 

00 

- . 1491689 iE 

00 

-.146890916 

00 

- , 1446446ie 

00 

-.142429656 

00 

-.140245696 

00 


.306173926 01 
•31056328E 01 
.315021786 01 
.319550616 01 
.324150986 01 
.328824146 01 
.333571346 Ql 
.338393876 01 
•34329304E 01 
.348270206 Ol 
.353326716 01 
•35846396E 01 
.363683386 01 
,368986426 ol 
.374374566 Ql 
.379849316 Ol 
,385412216 01 
.391064826 01 
,396808766 01 
.402645656 01 
,408577166 Ql 
.414604986 01 
•42073084E ol 
,426956526 Ol 
.433283816 01 
.439714546 Ql 
.446250576 Oi 
.452893826 Ql 
,459646226 Ql 
.466509756 ol 
.473486426 Ql 
.480578286 OX 
.487787426 Ql 
,495115966 01 
.502566076 Ql 
.510139956 01 
.517339846 Ql 
.525668026 Ol 
,533626826 Ql 
.541718596 Ol 
.549945726 01 
.558310676 Ol 
.566815916 01 
,575463966 Ol 
.584257376 01 
.593198756 Ol 
.602290736 Ol 
.611535996 Ol 
,620937246 Ol 
.630497246 ol 
.640218776 Ql 
,650104676 01 
.660157796 01 
.670381046 Ol 
•68077734E Ql 
.691349666 Ol 
,702101006 Ol 
.713034376 Ol 
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.lOOOOOOOE 01 
.lOOOOOOOfc 01 
.lOOOOOOOE 01 
•lOODOOQOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
•lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
•lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
•lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
,100000006 01 
•lOOOOOOOE 01 
•lOOOOOOOE 01 
•lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
•lOOOOOOOE 01 
.lOOOOOOOE 01 
,100000006 01 
.lOOOOOOOE 01 
•lOOOOOOOE 01 
.lOOOOOOOE 01 
•lOOOOOOOE 01 
.lOOOOOOOE 01 
•lOOOOOOOE 01 
•lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
•lOOOOOOOE 01 
.lOOOOOOOE 01 
•lOOOOOOOE 01 
.lOOOOOOOE OX 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
•lOOOOOOOE 01 
•lOOOOOOOE 01 
• lOOOOOOOE 01, 
•lOOOOOOOE 01 
•lOOOOOOOE 01 
.lOOOOOOOE 01 


.2080000QE 01 
.209000006 01 
.210Q0000E 01 
.21100000E 01 
.212000006 01 
.21300000E 01 
.2140D000E 01 
.21500000E 01 
.21600000E 01 
.217000006 01 
•21800000E 01 
.219QOOOOE 01 
.22000000E 01 
.22100000E 01 
'.22200000c 01 
•2230000UE 01 
.22400000E 01 
,225>0000QE 01 
.22600Q00E 01 
•22700000E 01 
.228-ODOO'OE 01 
.22900000E 01 
•23000000E 01 
.23100000E 01 
.23200000E 01 
,233000006 01 
,23400G00E 01 
.23500000E 01 
.23600000E 01 
.23700000E 01 
.23800000E 01 
.23900000E 01 
.24000000E 01 
.241000Q0E 01 
.2420000UE 01 
.243000006 01 
.24400000E 01 
.24500000E 01 
,246000006 01 
.24700000E 01 
•24800000E 01 
.24900000E 01 
.25000000E 01 
.2i5lUOOOaE 01 
.25200000E 01 
.2ij30a000E 01 
.264000006 01 
.256000006 01 
.2560000QE 01 
.257000006 01 
.25800000E 01 
.259000006 01 
.2600000UE 01 
,261000006 01 
•26200000E 01 
.2630QO00E 01 
.264Q0000E 01 
.26500000E 01 


-.13809239E 00 
-.135969426 00 
-.13387644E 00 
-.13181314E 00 
-.12977919E 00 
“.12777428E 00 
-.12579810E 00 
-.12385035E 00 
-.12193Q72E 00 
-.12003894E 00 
-.1181747UE 00 
-.11633772E 00 
-.11452772E 00 
-, 11,2744446 00 
-.11098756E 00 
-.109256916 00 
-.107852126 00 
1058729-76 00 
-.10421926E 00 
-.10259068E 00 
-.100987006 00 
-.994079886-01 
-.978533996-01 
-.963230046-01 
-.948165786-01 
-.933338876-01 
-.918747216-01 
-.904388596-01 
- . 89026090E-01 
-.87636203E-01 
-.86268998E-01 
-.849242636-01 
-.83601811E-01 
-.823014436-01 
- . 3102297QE-01 
797662046-01 
- . 785309616-01 
-.773170636-01 
7&124332E-01 
-.749525956-01 
7,i8016SlE-01 
-.726/14286-01 
-.715616646-01 
-.7 ;j4 72236E-01 
-.694029856-01 
-.683537576-01 
- .673244026-01 
-.66314772E-01 
-.653247236-01 
-.64354U3E-01 
-.6340280 16-01 
-.624/06596-01 
-,6155754''E-01 
-.606633386-01 
- . 59787905E-01 
- . 589311246-01 
- , 58092874E-01 
-.5;2730366-QI 


.724152846 Oil 
.735459496 01 
,746957406 01 
.758649726 01 
.770539586 01 
.782630156 01 
.794924596 01 
.807426096 01 
.820137836 01 
.833063026 01 
.846204826 01 
.859566436 01 
.873151026 01 
.886961736 01 
,901001706 01 
,915274026 01 
,929781766 01 
.944527936 01 
.959515516 01 
.974747396 01 
,990226436 01 
.100595546 02 
.102193696 02 
.103817366 02 
.105466796 o2 
.107142226 02 
.108843876 02 
.110571946 02 
.112326626 02 
.114108096 02 
.11591650^ 02 
,117751986 02 
.119614636 02 
.121504556 02 
.123421796 Q2 
.125366386 02 
.127338316 q2 
.129337566 02 
.131364046 02 
.133417666 02 
.135498276 02 
.137605676 02 
,139739626 02 
.141899856 02 
.144086026 02 
.146297746 02 
.148534556 02 
.150/95966 02 
.153081406 02 
,155390226 02 
.157721736 02 
.160075156 02 
.162449626 02 
.164844216 02 
.167257916 02 
.169689626 02 
.172138156 02 
.1/4602236 02 
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.iOOOOOODE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
-lOOOOOOOE 01 
-lOOOOOOOE 01 
-lOOOOOOOE 01 
-lOOOOOOOE 01 
-lOOOOOOOE 01 
-lOOOOOOOE 01 
-lOOOOOOOE 01 
.lOOOOOOOE 01 
-lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.100000006 01 
.100000006 01 
.lOOOOOOOE 01 
•lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
•lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
,100000006 01. 
•lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
.lOOOOOOOE 01 
-lOOOOOOOE 01 
-lOOOOOOOE 01 
-lOOOOOOOE 01 
.lOOOOOOOE 01 


.26600000E 01 
.26700000E 01 
.26800QOOE 01 
.26900000E 01 
-27000000E 01 
-271C0000E 01 
,272000006 01 
.273000006 01 
-27400000E 01 
.27500000E 01 
.27600000E 01 
.27700000E 01 
.27800000E 01 
.27900000E 01 
,260000006 01 
.28100000E 01 
,282000006 01 
,283000006 01 
.284000006 01 
-285000006 01 
.286GOOOOE 01 
-2870000UE 01 
.288000006 01 
,289030006 01 
.290000006 01 
,291000006 01 
.29200000E 01 
.293000006 01 
.294000006 01 
.295000006 01 
.29600000E 01 
,297000006 01 
.298000006 01 
,299000006 01 
.300000006 01 
.301000006 01 
.302000006 01 
.303000006 01 
.304000006 01 
.305000006 01 
,306000006 01 
.307000006 01 
.308000006 01 
.'309000006 01 
.310000006 01 
.311000006 01 
.312000006 01 
.313000006 01 
.314000006 01 
,315000006 01 
case C0MPLETED 


", 564/14966^01 
'•.556881 4 06 "01 
'•,549228576-01 
541/55426-01 
- . 5.1446088E-01 
'..52734394E-0X 
“-52040359E-01 
- .51363888E-01 
-.50704886E''01 
- , 50 06326 0E»01 
--49438922E-01 
-.48831784E-01 
“.46241763E-01 
■“-4/668776E-01 
471127456-01 
465/3592E-01 
"-46051242E-01 
-. 45545623E-01 
». 4505666/6-01 
-,44584304£-0l 
-,441284716-01 
-,43689103E«01 
•-43266142E-01 
“-42859529E-01 
-.4246920/E-01 
-.42095124E-01 
-.417372286-01 
--41395470E-01 
-.410698036-01 
- . 407601836-01 
404665686-01 
40188916E-01 
--39927191E-01 
--39681356E-01 
-.394513796-01 
-,3923722/E-OI 
-.390388716-01 
- . 38856285E-01 
-,38669444E-0X 
-.385383246-01 
~,38402906E-01 
- . 38283171)6-01 
-.381791016-01 
--38090685E-01 
- . 38Q17908E-01 
-.3/960/62E-Q1 
-.3/9192396-01 
-.37893333E-01 
-.378830416-01 
- . 3/S8836UE-01 


.177080496 02 
.17957145E 02 
-18207356E 02 
.184585146 02 
.187104436 02 
.18962956E 02 
-19215855E 02 
.194689316 02 
.197219656 02 
.199/47286 02 
•20226978E 02 
.204/8465E 02 
.207289276 02 
.20978Q92E 02 
.212256796 02 
.214713956 02 
.21714941E 02 
.21956006E 02 
.221942746 02 
.224H9418E 02 
.226611076 02 
.22889003E 02 
.231127616 02 
.233320346 02 
,235464726 02 
,237557206 02 
,239594256 02 
.241572336 02 
.243487906 02 
-24533746E o2 
.2471175/6 02 
.248824826 02 
.25045589E 02 
,25200/516.02 
.253476576 02 
-25486001E 02 
.25615495E 02 
.257358626 02 
•25846843E 02 
.259481966 02 
.260396966 02 
.26121139E 02 
.2619234QE 02 
.262531386 02 
.263033936 02 
,263429906 02 
.263718376 02 
.26389866E 02 
.a6397036£ 02 
,263933306 02 
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ANUM = AHI*DI + (l.0-B**2)*AN2*D2 
DEH = DI*DI+{I.O - B*»2)*D2»D2 
RE = ANUM/ DEN 
REC = -I .0/N3 



FORM Kl EQUATION SET 

DO 380 11- I, IQ 
DO 380 IJ = I, IK 
RK = I.O/AK(IJ) 

DBN = AN1*DI +BSQ»AN2*D2 

Q8 = 1.0 - Q (II)*BM1 

BQ = Q(l l)*BSQ 

DN = ANI*D2 - AN2»D1 

DB = DI»DI + BSQ»D2*D2 

TOP = DBH*QB + BQ*DN + RK»DB 

AQB = Q( I l)»BMI 

BN = ANI*ANI + BSQ*AN2*AN2 

BOT = AQB* BN 

AKP = T0P/80T 






FORM QUADRATIC EQUATION IN K 
00 399 II = 1, IQ 

AA = Q(l l)*BMI*(ANI*ANI + BSQ»AN2*AN2) 

BB = (ANI*DI + BSQ«AN2*02)*(I.0 - Q(II)*BMI 
CC = (DI*0I + BSQ*D2*02) 

DISC = SQRT (BB*BB - i|.0*AA*CC) 

AKI = (-BB + 0ISC)/(2.0*AA) 

AK2 = (-BB - DISC)/(2.0*AA) 
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